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Abstract 

The reliable simulation of extensive air showers induced by different primary particles (e. g. 
proton, iron, gamma etc.) is of great importance in high energy cosmic ray research. The 
CORSIKA pQ is a standard Monte-Carlo simulation package to simulate the four dimensional 
evolution of Extensive Air Shower (EAS) in the atmosphere initiated by gamma, hadrons 
and nuclei. CORSIKA has different high energy interaction models like DPMJET, QGSJET, 
NEXUS, SIBYLL, VENUS and EPOS which are based on different theoretical frameworks. 
The influence of different hadronic interaction models, viz., QGSJET and DPMJET on the 
lateral distribution functions (LDF) and muon to electron ratio of cosmic ray EAS induced 
by 10 17 eUto 10 20 eV proton and iron primaries are explored in this work. 
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1. Introduction 



The flux of cosmic ray particles extends to energies far beyond the reach of any Earth-based 
particle accelerator. Again direct observation of the first interaction of UHE cosmic rays in the 
upper atmosphere is impossible due to very low flux. The information on the first interactions 
can be derived from the obtained characteristics of primary induced cascades of secondary 
particles, known as extensive air showers(EAS); which requires ground based large detector 
setup. So the interpretation of EAS characteristics for primary energy above the energy at- 
tained by man made accelerator is dominated by model predictions and detailed simulations. 
Differences and uncertainitties in these hadronic interaction models play the crucial role in 
the prediction of EAS characteristics. For air shower interpretation the understanding of mul- 
tiparticle production in hadronic interactions with small momentum transfer is essential [2]. 
Due to energy dependence of the strong coupling constant as, soft interactions cannot be cal- 
culated within QCD using perturbation theory. Instead, phenomenological approaches have 
been introduced in different models incorporated in CORSIKA. Study of LDF(Lateral Distri- 
bution Function) for muons and electrons/positrons and photons is important for low energy 
hadronic model dependence study, especially LDF of muons. [5] . Again shower maximum 
prediction varies with defferent interaction models. Models that predict a higher scattering 
cross-section, higher multiplicity produces shower maximum at smaller slant depth. In the 
energy range 10 17 eFto 10 20 eFlots of ambiguities regarding the composition of primary cosmic 
rays are prevailed. Therefore simulations for extreme assumptions for primary (viz. protons 
and iron nuclei) are done. Corresponding predictions for different observables are calculated 
for different interaction models. The measured data should be in between the predictions 
obtained for the extreme assumptions. If the data are outside the proton-iron range for an ob- 
servable, there is is an indication for an incompatibility of the particular hadronic interaction 
model with the observed values, jl]. 

Hadronic Interaction Model: 

In order to explain hadron-hadron/nucleus- nucleus collision above the attainable energy in 
accelerator one has to heavily rely on proposed hadronic interaction models. QGSJET 
(Quark Gluon String with JETs): [3] is the most successful model in describing the high 
energy hadronic interactions. This model offered relatively easy approach to the simulation 
of cosmic ray interactions at higher energies, as well as ensured a good agreement to the 
accelerator data at lower energies [6]. This QGS based model treats the hadronic and nuclear 
collisions in the framework of Gribovs reggeon approach as multiple scattering processes, where 
individual scattering contributions are described phenomenologically as Pomeron exchanges. 
The Pomeron corresponds to microscopic partons (quark and gluon) cascades, which mediate 
the interaction between the projectile and the target hadrons, and consists of two parts: soft 
and semihard Pomerons. So in this model the appropriate cross sections for the inelastic 
interaction between hadrons i and j are calculated using the expression [1], 



where Xij( s >b) and Xij( s >b) represent the soft and semihard pomerons respectively, s is 
the centre of mass energy and b is the impact parameter. 
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The total cross section in this model is given by the formula [2], 

= ^ / d% [l - exp 6) + &))}] [2] 

where is so-called the shower enhancement co-efficient, for pp interactions epp = 1.5. 

DPMJET (Dual Parton Model with JETs): is based on the two components Dual 
Parton Model [7] and contains multiple soft chains as well as multiple minijets. As QGSJET, 
it relies on the Gribov-Regge theory (8J |9] and the interaction is described by multi-Pomeron 
exchange. Soft processes are described by a supercritical Pomeron, while for hard processes 
additionally hard Pomerons are introduced. High mass diffractive events are described by 
triple Pomerons and Pomeron loop graphs, while low mass diffractive events are modelled 
outside the Gribov-Regge formalism. The total cross section in this model is given by [3], 



a t ij (b,s)=4irCbdb[l-exp( X i j (sM [3] 
where Xij{ s ify represents total Pomeron in this model. 



Method 



Here we have used CORSIKA-6735 Simulation code to generate EAS using two hadronic 
interaction models viz. QGSJET-01 and DPMJET. We have considered two primary masses 
(proton and iron) and four primary energies(10 17 e V,10 18 e V,10 19 eF and 10 20 eV). There are 
altogether 16 sets of events with 1000 showers each. Simulations are performed for vertical 
showers. GHEISHA sub routine and THIN options are selected for flat horizontal detector. 
Data are first extracted by the inbuilt FORTRAN program and then these are analyzed using 
C++ program in the ROOT environment. 




Results: 

In order to study model dependence, we have chosen the following parameters: 
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(i) the ratio of total number of muons to primary energy as a function of primary energy 
are plotted in Fig.l The predictions by the DPMJET model for both the primaries have higher 
values except for proton induced shower at 10 17 eV. Slopes of the curves for DPMJET and 
QGSJET are almost equal. 

(ii) Average depth of shower maximum < X max > is plotted as a function of logarithm of 
primary energy as shown in Fig. 2. The secondary data from Hires experimental setup are also 
superimposed into it. It is observed that irrespective of High Energy Hadronic Interaction 
Models considered, data are in between the predictions by proton induced showers and the 
same by iron induced showers. 

(iii) Ratio of muon no to electron number as functions of primary energy and shower size 
are shown in Fig.3 & Fig.4. Both these figures show that muon number is more for iron 
induced showers as compared to that produced for proton induced showers. However there is 
no significant difference so far as the two models are concerned. 

(iv) Lastly, the density of muons and electrons at different core distance are calculated from 
the simulated data in 2m step and their ratio is plotted as function of core distance for proton 
and iron primaries. Fig. 5 shows these results for two hadronic interaction models considered 
here. It is seen that there are more fluctuations at larger distances and no significant model 
dependence in all the cases. 




Conclusion 

From the study of the different parameters it is found that the ratio of muon number to 
primary energy is best suited for distinguishing the effect of interaction model. DPMJET is 
found to produce more muons per GeV primary enegy compared to QGSJET-01 model. 
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Figure 5: Lateral distribution of muon to electron ratio for proton and iron primaries 
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